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Abstract 

We examine cosmological constraints on the lepton number breaking scale 
in supersymmetric singlet majoron models. Special attention is drawn to the 
model dependence arising from the particular choice of a certain majoron 
extension and a cosmological scenario. We find that the bounds on the sym- 
metry breaking scale can vary substantially. Large values of this scale can 
be allowed if the decoupling temperature of smajoron and majorino exceeds 
the reheating temperature of inflation. In the opposite case an upper bound 
depending on the majoron model can be obtained which, however, is unlikely 
to be much larger than 10 10 GeV. 



Singlet majoron models Q provide the simplest extension of the standard model 
with spontaneously broken lepton number and neutrino masses generated via the 
celebrated see-saw mechanism. Therefore these models are of particular interest 
for cosmology. The well-known cosmolgical constraint on the conventional singlet 
majoron model comes from the late-decay of heavy neutrinos into majorons (J) and 
lighter neutrinos. In the model with one singlet S the amplitude for this decay has a 
strong suppression ||. Then the standard cosmological constraint derived from the 
closure density of the universe |§ bounds the right handed neutrino (RHN) mass to 
be < 10 3 GeV taking the largest possible neutrino mass of 30 MeV and a mixing 
angle of order 10 -4 . For generic models with an arbitrary number of singlets the 
bound can be shifted to ~ 10 9 GeV. Clearly, for decreasing neutrino masses these 
bounds become weaker and for masses below ~ 40 eV no information about the 
RHN scale can be obtained at all. 

In supersymmetric models, cosmological considerations provide independent in- 
formations due to the presence of the superpartners of the majoron (the smajoron a 
and the majorino if)) as discussed first by Mohapatra and Zhang |4j]. Because their 
coupling is surpressed by the lepton number symmetry breaking scale Vj, they can 
be expected to decouple when relativistic. In comparison to supersymmetric axion 
models || the partners of the majoron tend to couple much weaker to the standard 
particles due to the lack of an anomaly and the smallness of the neutrino masses. 
One might therefore ask if the constraints on Vl from baryogenesis and overclosure 
will be more restrictive than corresponding constraints on axion models. We con- 
sider the range from ~ 10 10 GeV to 10 12 GeV as particulary attractive for Vl since 
several independent hints from the MSW solution of the solar neutrino problem, the 
hidden sector SUSY breaking and axion physics point to such scales. A question 
of particular interest might therefore be to what extent such large values of Vl are 
compatible with cosmological restrictions on SUSY majoron models. 

In the below we examine this matter taking two kinds of model dependences into 
account. On the one hand we discuss the freedom in the field theoretical parameters 
relevant for cosmology like masses and couplings of the smajoron and the majorino. 
On the other hand we analyze the influence of the cosmological scenario, namely the 
influence of inflation which for large Vl may be able to wash out the relic densities 
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of smajoron and majorino. 



We start by briefly introducing the class of models we are considering here. Their 
superpotentials which come in addition to that of the MSSM can be written as 

W = h ij L i HN j + f ija NiNjS a + W(S a ) (1) 

with three families iV; of RHNs and gauge singlet fields S a . For most of our discus- 
sions soft breaking terms of the standard type 

Vsoft = mj /2 £ |0 a | 2 + m 3/2 <Pa^~ + (A- 3)W + h. a) (2) 

are assumed. We have in mind that supergravity is broken in the hidden sector. 
The whole potential is invariant under Ul(1) with Ql{L) = 1, Ql(N) = —1 and 
Ql(So) — la, which is spontaneously broken by the VEVs v a of the gauge singlet 
fields S a . Apart from being larger than m 3 / 2 the Ul{1) symmetry breaking scale 
V L = (J2 \ v a\ 2 QlY^ 2 is taken as a free parameter. The freedom for model building 
consists in the number of singlet fields S a and the choice of their potential W(S a ). 
In this paper we assume that R-parity is unbroken. 

For cosmological considerations we extract now the model-dependence of masses 
and lifetimes of smajoron and majorino. Let us first examine the interactions of the 
majoron supermultiplet $ = ((a + iJ)/\^2,?p). The interaction with the standard 
particles is given by the usual form 

M f 

W = hLHN + —NN + J -NN<& (3) 

with the generation indices suppressed. We also need to know higher-dimensional 
effective interactions. The one most important for our purpose comes from a one 
loop diagram built up from the interactions LHN and <&NN. It leads to the effective 
D=5 operator ($ + <I)LL. Putting the fields on shell the diagram vanishes so that 
there is no decay mode of $ with a surpression. However, attaching a gauge 
field V gives the coupling 
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Here h denotes the dominant one of Yukawa couplings in eq. ([!]) . This diagram 
can keep the majoron supermultiplet in equilibrium after the phase transition. 

To describe the interactions among the components of $ which generally depend 
on the form of W(S a ), the nonlinear approach is easier to handle since the majoron 
and the majorino are completely rotated away from the superpotential as well as 
from the soft terms. Only the Kahler part needs to be considered. Starting with 
the minimal Kahlerian K = J2a S a S a and writing S a = v a exp(g ffi $/Vi,), we arrive at 
the following Lagrangian in component fields : 

\plx 11- 
L = (1 + -y—rfi-dpadpa + -d^Jd^J + itpa^tp) 

- — (d^Jtya^ip + (4 fermion terms) + 0(l/v 2 ) (5) 
with x = -TTT^a ■ (6) 

a * L 

The model dependence is now condensed in the value of x which influences the decay 
rates for o —>■ 2J and o — > 2ip. 

The value of x can be of order one or zero at tree level. Even in the second case 
a nonzero value can be generated by radiative corrections. To discuss this, let us 
take the simplest model. It is specified by the superpotential 

W = hijL ai H a Nj + fijNiNjS + X(SS' - /i 2 )F (7) 

which results in the following soft terms for the gauge singlet fields : 

V so ft= m Jl0| 2 + m 3/2 A n fnns + A s Xss'y - A y X/j, 2 y + h. c. . (8) 

<f>=n,s,s' 

The boundary condition to be fulfilled at Mpi is = TO3/2, A n = A s = A and 
A y = A — 2. The value of x for the above potential is x = (v 2 — v' 2 )/(v 2 + v' 2 ). 
Minimization leads to 

m 2 , — m 2 m 2 , — m 2 , nS 

r, = ~ — (9) 

2\ 2 y 2 + m 2 s + m 2 , 4m 2 3/2 ' 1 ; 

with y = (A y — A s )m^/2/2\. Obviously degenerate soft masses result in x = which 
is not surprising since in that case a permutation symmetry s <-> s' exists. The 
only way to break this symmetry is to take m s 7^ m s i which can be achieved in two 
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ways : One can assume nonuniveral mass terms at Mpi thereby departuring from the 
standard soft terms (Q) or one can consider renormalization group effects. The first 
option clearly allows to produce any value of x. So let us discuss the implications 
of RGE starting with universal boundary masses. Due to the asymmetric coupling 
of s and s' to the RHNs their masses renormalize differently and one finds || 



f 2 „ , ^ (M. 



x ~ -^—(1 + A z )ln 



pi 



64tt 2V ' V V L 

.os I 9 1 /l0 12 GeV\\ , , 

= />(l + ^)(2.2xlO- + ^>n^^JJ . (10) 

So it appears to be difficult to get x > 0.1. The lifetime of the smajoron decay 
a — ► 2 J is given by 

= / V L \ /l0 2 GeV\ 3 , s 



In another model discussed in ref. [|7| only one change appears. The value of y 
is negligible, and so the value of x roughly doubles. 

The majoron has mass zero forgetting about possible ?7L(l)-breaking gravity 
effects. Considering this effect appearing in the D=5 operators it was found that 
the cosmological mass density constraint gives a strong bound on the lepton number 
breaking scale : Vl < 10 TeV ||. In this paper we neglect this effect. The smajoron 
cannot escape receiving a mass m a ~ 0{m$/2) from soft terms. The majorino tree 
level mass is model-dependent and can take values m^ ttree ~ 0((m 3 /2/VL) k m 3 /2) 
with k = or 1 typically. A small tree level majorino mass (k > 0) is e. g. realized 
in the model of ref. [7|][] or in the above model for nonstandard soft terms satisfying 
A s = A y at Mpi. However, in analogy to the heavy quark axion ||, it receives a 
model-independent one-loop correction m^j oop = (l/167r 2 )f 2 Am 3 /2, where / is the 
dominant one of the couplings fij a . Additionally, there can arise another radiative 
mass due to the RG evolution of the trilinear soft couplings in such models. The 
induced mass depends on / and the completely unrestricted couplings in the W(S a )- 
sector [[IT]] and is therefore quit model-dependent. In any case the majorino mass 
lr This model can provide a light majorino with the standard universal soft terms. But it has 



another zero mode than the majoron, which may cause additional cosmological impacts 10 
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can take values far below m.3/2. To illustrate the cosmological implications we will 
use the above model-independent value m^^ oop and examine the cosmological effects 
depending on the value of /. 

When the temperature of the universe falls below the scale Vl, smajoron and 
majorino can be kept in equilibrium via effective interactions with light fields. Con- 
sidering the D=5 effective interaction in eq. smajoron and majorino get out of 
equilibrium below the temperature 

T L = 6.2 x 10 17 /fVW 100 (— V GeV . (12) 

\ Trip j 

Hereafter we introduce the notation f n = 10 n f. The actual decoupling temperature 
is the lower one between Vl and Tl- From this we conclude that smajoron and 
majorino decouple when they are relativistic. Such hot relics can cause the problem 
of the excessive energy density if the particles are late-decaying. In regards of infla- 
tion this problem should be examined in two ways. The usual see-saw mechanism 
assumes a large scale Vl ~ 10 12 GeV for the lepton number breaking. It is bigger 
than the reheating temperature Tr which should be less than about 10 10 GeV in 
order to avoid the gravitino problem | H| . In such a case (Vl and Tl > Tr) the 



relic density of majorino or smajoron can be washed away. However, there can be a 
significant amount of regerated relics after the reheating. We will later consider the 
effects of regeneration. 

It can also happen that smajoron and majorino decouple after inflation (Vl 
or Tl < Tfj). For this case we note that the crucial model-dependence of the 
cosmological problem resides, especially, on the values of x and m^. We begin with 
the discussion of this case. The ratio of the relic density to the entropy density 
is then of order 10~ 3 . In the case that the smajoron decays into two majorons, 
the decay-produced majorons should be red-shifted away sufficiently in order not to 
upset the standard nucleosynthesis prediction. The majorons from smajoron decay 
should not supply an energy density larger than three tenths of the density of one 
neutrino species, which implies 



nin \ 1 ~ n \ 1/2 . _ _ 5 



10 2 GeVy Vlsec 



<8.1xl(T 5 . (13) 
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From the lifetime given in eq. ( IT[) one finds 



V L < 1.8 x 10 w x 



(14) 



Thinking of models with nonuniversal scalar masses one may expect x = 0(1) so 
that it seems difficult to push Vl higher than 10 10 GeV. Taking the universal bound- 
ary conditions for the soft terms the value of x depends on the Yukawa coupling 
/ between majoron superfield and RHNs. For the two explicit models we have 
mentioned here the situation can be summerized as 



When the value of / becomes extremely small, the decay channel of smajoron into 
two majorons is suppressed and those into two neutrinos or two sneutrinos become 
important. This is the region of 



where the constraint in ref. |4[] is applicable. 

A more severe constraint comes from the decay of majorinos. Being R-parity 
odd, each majorino will eventually produce at least one LSP (ordinary neutralino 
X°) whose mass in generally exceeds the GeV-range. Then those relics of secondary 
LSP cause an overclosure of the universe. To avoid this problem the majorino should 
decay before the decoupling time of neutralinos which is about 10 _6 (20 GeV/m x o) 2 
sec. If the majorino is heavier than the sneutrino the decay ip — > vv produces a 
strong bound : 



One can also find a bound of a similiar order of magnitude if the decay channel 
ip -> uuH° |§ is allowed. 

Of course the majorino can be the LSP istself in which case it is a warm dark 
matter candidate. Such a majorino has to fulfill the overclosure bound < 2 
keV which translates to a restriction on / 




(15) 





(16) 




(17) 
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The model has now to be such that the tree level mass is below ~ 2 keV. Assuming 
that x comes entirely from RG effects eqs. (17, 15) give 



which is quite restrictive. It should, however, be stressed that this bound is influ- 
enced by the magnitude of RG effects and is therefore model dependent. Taking a 
model with a light majorino and nonuniversal soft masses it would be completely 
invalidated and provided / satisfies eq. fllTD the strongest constraint on V L would 
be given by eq. flU]). 

Let us now consider the case that smajoron and majorino decouple before or 
during inflation (Vj, and Tl > T R ). As already mentioned our concern is the regen- 
eration of relics. In terms of the ratio of relic density to entropy density Yx = nx / s 
the Boltzmann equation of the regenerated population for a particle X becomes 

dY x (Z T v)n 2 r 



dT sHT 



(19) 



where we neglected the decay terms. We can explicitly integrate the Boltzmann 
equation from the reheating temperature T R to some temperature T with an initial 
condition Y(T R ) = 0. If T R > T, we get 

Yx(T) * ( ^^™ 2 = 5.9 x 1CT 6 (E T v)M Pl T R . (20) 
H(T R )s(T R ) 

We see that Yx is determined by (Htv) and T R . For our case, the D=5 interaction 
in eq. (§]) is most important for the regeneration and (Etv) ~ g 2 h A /16 3 ir 5 Vl. Then 
we obtain the regenerated population for smajoron or majorino : 

Y a ^ 6.5 x IO-VKt^tV (t^) ■ (2D 



VlOeV/ V10 10 GeV 

We first consider the relevant constraints coming from the regeneration of ma- 
jorinos. If the majorino is stable or its lifetime is longer than the age of the 
universe t , the relation 
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has to be fulfilled in order not to overdose the universe. This puts a bound on the 
neutrino mass 

m - <17/ ' HH {—) eV - (23) 

If the majorino is unstable and < t , the overclosure density bound is applied to 
the decay-produced neutralinos. For this we just replace by m x o in the above 
equation. 

There is another contstraint coming from the high-energy neutrinos when the 
majorino (smajoron) is unstable. The decay of smajoron and majorino generally 
produces high-energy neutrinos which can be observed by experiments. For the 
energy range of our interest, the best bounds come from the experimental upper 



limit on upward-going muons from the 1MB detector ||12|| . Here we quote the rough 
order of estimation 



2 



18 


r to 




W 


18 


( Tip 




Kt 



4/3 



< 10" 18 m T t >t a . (24) 

Considering the decay of the majorino into a neutrino and a sneutrino, we find for 

V L < 1.6 x f^V' 4 f!5!^XV /4 GeV, (25) 



VlOeVy \ rriip J \ T, 



R 



and for > to 

Here the majorino lifetime is given by 

,n tniB / Vl \ 2 /10eV\ 2 A0 2 GeV\ 

The above results are summarized in the figures. Fig. 1 shows the forbidden 
region in the Vl — space when the decay ip — > vv is allowed. The lines (ii) or 
(iii) are very sensitive to the value of / and they provide no constraint unless / is 
larger than about 10~ 4 . As one expects, a rather strong constraint can be found for 
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larger values of / and Tr. We note that Vl is limited between Tr and the value in 
eq. fl25|) in order to have, e.g., a 40 eV neutrino as a dark matter candidate. If the 
reheating temperature is pushed up to 10 10 GeV, the line (iv) goes down below the 
line (a) and we find a strong bound on the neutrino mass: 

,„,-2«/10 2 GeVV /2 /10 1 °GeV\ 1 '" i „ 

m '< i8 ''"(— ) HH eV - (28) 

The size of the constrained region becomes smaller if the lifetime of the majorino 
is fairly long or if it is stable. When majorino is stable, only the lines (ii) and (iii) are 
to be taken. As an example for a long-lived majorino which can decay into the LSP, 
let us take the decay ip - > vvH Q with the Higgsino H° as the LSP. The lifetime can 
be found in ref. 0] and together with the constraint from the high-energy neutrinos 
it results in the lower part of Fig. 2. The constrained region becomes smaller for 
larger Tr as the curve moves below the line (a). For instance, there is no region 
forbidden if Tr = 10 10 GeV. 

Regenerated smajorons play a less important role than majorinos. With x ~ 1 
one finds no constraints since then smajorons decay fast enough. In a model with the 
RG-induced x ~ 10 -2 / 2 , the smajoron decay into two neutrinos can be important 
for small values / < 10 -3 , for which the constraint from the high-energy neutrino is 



the same as in eqs. (25, Eq). The effect increases for smaller values of / and larger 



values of Tr. This behavior is shown in the upper curve in Fig. 2. 

In conclusion we have discussed the cosmological constraints on SUSY majoron 
models. The majoron superpartners decouple when relativistic and we should dis- 
tinguish two cases, the decoupling temperature is lower or higher than the reheating 
temperature after inflation. 

In the first case we identified the mass of the majorino and the quantity x 
which determines the strength x/Vl of the decay o — > 2 J as the crucial parame- 
ters influencing the cosmological constraints. From the decay of the smajoron we 
obtained a bound Vl < xlO 10 GeV where in models with large RG effects or nonuni- 
versal scalar soft masses one can think of x as being of 0(1). Assuming standard 
soft scalar masses, however, in our examples x was determined by small RG effects 
leading to typical values x ~ 10~ 2 / 2 where / is the coupling of the majoron to the 
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RHNs. Stronger bounds were found from the majorino decay. If the majorino can 
decay into the LSP, Vl is constrained to be < O(10 2 GeV). However, there exist 
models with a negligible tree level majorino mass and a radiative mass contribu- 
tion proportional to /. In such models the majorino can be the LSP and to avoid 
overclosure / has to fulfill / < O(10 -3 ). If at the same time x is generated by RG 
effects, a small / surpresses the decay a — > 2 J which results in a typical bound 
Vl < O(10 3 GeV). It should be stressed that this bound is avoidable : A model with 
a keV majorino and x of order one would be fully consistent up to Vl ~ 10 10 GeV. 

For the second case we examined the effects of the regenerated relics. The con- 
strained region in the Vl — m u plane is strongly dependent on the values of / and 
the reheating temperature Tr as well as the decay channel of the majorino. So one 
needs more informations on the various parameters to draw some definite conclu- 
sions. Capitalizing the general pattern, we find a relatively strong constraints if the 
majorino decays into a sneutrino and a neutrino is allowed. In this case a higher 
value of Vl together with a 10 eV neutrino can be allowed for a small value of the 
Yukawa coupling / < O(10 -3 ). If the above decay channel is not allowed, we con- 
clude that a higher value of Vl can be easily achieved. In either case the bounds on 
heavy neutrinos tend to be more restrictive. We finish by making a remark on the 



recent work |R| which showed that the supersymetric singlet majoron model can 
explain both chaotic inflation and baryogenesis. In this scenario the inflaton is iden- 
tified with the scalar partner of a RHN. Then the primordial density perturbation 
fixes the RHN mass at M ~ 10 13 GeV. Furthermore the reheating temperature of 
order 10 10 GeV could be obtained with reasonable values of parameters. Therefore 
this scenario should be such that the decoupling of majorino and smajoron occurs 
before the reheating and should follow the constraints we have investigated. For 
instance, assuming a Yukawa coupling / of order one, the above discussion shows 
that the majorino should be lighter than the sneutrino. 
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ure captions 

Fig. 1. For the decay ip — > vv. The curves are normalized to T R = 10 8 GeV 
and / = 0.1. 

- The lines (a) and (b) define the region where V L > T R and T L > T R , 
respectively. The region below the line (tl) corresponds to r(ip — > vv) < 
to- 

- Below the line (i) the yukawa coupling h remains in perturbative region 
h < 1, that is, 

(io*lev) (iTiv) <6 - 1//l - 

- The lines (ii) and (iii) arise from the overclosure bound for (quasi-) stable 
majorinos and decay-produced neutralinos, respectively. The region left 
to the lines is allowed, (cf. eq. (23) in the text.) 

- The lines (iv) and (v) show the constraints from the decay-produced 
energetic neutrinos. The region below the line (iv) and above the line (v) 
is allowed, (cf. eqs. (25, 26)) 

Fig. 2. For the decay ip — > vvH° and a — > J J, vv. The curves are normalized 
to T R = 10 8 GeV and / = 10" 4 . 

- The lines (a), (b) and (i)-(v) are as explained in Fig. I. The lines (ii) 
and (iii) are now outside the experimental limit of m u . 

~ The lines (si, s2) represent r(a — > JJ, vv) = to- The line (t2) corre- 
sponds to r{jp — > vvH°) = t . 

- The lines (vi) and (vii) show the constraints from the decay-produced 
energetic neutrinos as in eqs. (25, 26). The region below the line (vi) is 
defined by 

„ /19 in8 ,-7/8 / m u y/ 8 ( y/ 8 ^ 10 8 GeV V /16 

and the region above the line (vii) by 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403217vl 



This figure "figl-2.png" is available in "png" format from: 
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